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■In  this  paper  we  develop5Jan  algebraic  approach  to  the  multiple  time  scale 
analysis  of  perturbed  linear  systems  based  on  the  examination  of  the/ Smith 
form  of  the  system  matrix  viewed  as  a  matrix  over  a  ring  of  functions  in  the 
perturbation  parameter.  This  perspective  allows  us  to  obtain  a  strengthened 
version  of  the  results  of  £3r}- and  to  provide  a  bridge  between  these  complex 
but  general  results  and  previous  explicit,  concepu tally  simple,  but  somewhat 
restrictive  results. such  as  described  in  [1],  [2].  ^In  addition,  -our  algebraic 
framework  allows -tfcs  to  investigate  a  variety  of  other  problems.  In  this  paper 
a#er  study  the  problem  of  developing  valid  time  scale  decompositions  in  cases  in 

/  which  weak  damping  .terms  discarded  in  the  approaches  in  {4} - [3}  must  BeT 

retained.  Also,  our  approach  exposes  the  role  of  the  invariant  factors  of  the 
system  matrix  in  determining  its  time  scales.  This  leads  naturally  to  the 


■  I-  s 


problem  of  time  scale  modification,  i.e.  invariant  factor  placement,  via  state 
feedback.^'  We -present  a  result  along  these  lines,  t  t 
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I.  Introduction 


This  paper  is  concerned  with  the  multiple  time  scale  analysis  of  the 
perturbed  N-dimensional  linear  system 

x(t)  =  A(e)x( t)  (1.1) 

where  A(e)  has  a  Taylor  expansion  in  the  small  parameter  a.  If  there  is  a 
drop  in  the  rank  of  A(e)  at  e  =  0,  the  system  (1.1)  is  termed  singularly 
per turbed  and  can  exhibit  multiple  time  scale  behavior.  The  analysis  of  such 
behavior  has  been  the  subject  of  a  number  of  previous  investigations.  In 
particular  several  researchers  [1],  [2],  [6],  [10],  [11],  [12]  have  made 
numerous  important  contributions  by  investigating  systems  in  what  we  will  call 
explicit  form^: 


w 

0* _ 

x2(t) 

= 

A11  A12  xl(t) 

fcA21  eA22  x2^ 


(1.2) 


fAll  A12l 


|A21  A22 I 


and  let  A 22  denote  the  Schur  complement  of  A^  in  A: 


(13) 


A  _  A  -  A  A  A 
22  "  22  *21  11  12 


(1.4) 


*The  forms  actually  considered  in  [1],  [2]  have  a  appearing  on  the  left-hand 
side  rather  than  the  right-hand  side.  There  is  no  significant  difference  in 
considering  the  form  (1.2),  since  the  systems  in  [1],  [2]  can  be  brought  to 
this  form  by  the  change  of  time  scale  t  =  t/a. 


It  is  known  that  if  A^  and  are  nonsingular,  the  eigenvalues  of  (1.2) 
occur  in  two  groups,  one  being  of  order  1  and  lying  "close"  to  the  eigenvalues 

iV 

of  A^,  and  the  other  being  of  order  e  and  close  to  the  eigenvalues  of  eA^ • 

If  both  the  latter  matrices  are  Hurwitz,  then  the  system  exhibits  well-behaved 

two-time-scale  structure,  in  the  following  sense: 

x  (t)l  fxlf(t)  +  x-  (et)  +  0(e)l  (1.5) 

=  1  iS  .  t  >  0. 

x2^ c)  x2s^6t^  +  °(e) 

where 

xif(t)  =  Aj  jX^(t)  ,  Xj£,(0)  =  x^(0)  +  Aj  ^  A12x2^) 

xls^>  =  "  A11  A12X2s^’  (16) 

and 

x2s^  =  A22  x2s(t)‘  X2s^°^  =  x2^0^- 


The  subcripts  s  and  f  denote  slow  and  fast  subsystems. 

The  0(e)  terms  in  (1.3)  are  uniform  in  t  >  0,  so  that  (1.5),  (1.6) 
provide  a  uniform  approximation  of  the  state  transition  matrix  of  (1.1).  That 


is. 


lim  sup  ||eA^e)t  -  T  le  ^  T | |  =  0 
elO  t>0 


(1.7) 


where 


The  decomposition  provided  in  (1.5)  -  (1.6)  or,  equivalently,  in  (1.7)  - 

(1.9)  has  found  significant  applications.  One  important  limitation  of  these 

results,  however,  is  the  assumption  that  the  system  is  given  in  the  explicit 

form  (1.2)  or  its  obvious  generalizations  (e.g.  by  expanding  the  A-matrix  in 

2 

(1.2)  to  include  a  third  row  of  blocks,  each  of  which  is  multiplied  by  e.  ) . 

On  the  other  hand,  there  are  several  advantages  if  the  system  has  the  form  in 

(1.2) .  Specifically,  there  is  a  simple  check  to  see  if  the  system  has  a  time 
scale  decomposition  in  the  sense  of  (1.7),  (1.8)  (namely  A. ^  and  A?9  must  both 
be  Hurwitz),  one  immediately  knows  what  the  time  scales  are,  and  the 
subsystems  describing  the  behavior  at  each  time  scale  are  easily  obtained. 

In  contrast  to  the  results  just  described,  we  have  the  work  of  Coderch, 
et  al.  [3]  which  had  as  its  principal  objective  the  development  of  a  general 
procedure  for  determining  if  the  system  (1.1)  has  well-defined  time  scale 
structure  and  for  constructing  a  decoupled  time-scale  decomposition  as  in 
(1.7)  with 

k  v  k 

Ad(e)  =  diag  (e  XAj  .  e  Ag . £  (1.10) 

(and  with  an  appropriate  choice  for  T)  wi thout  assuming  that  the  system  is  in 
the  special  form  of  (1.2).  This  objective  is  achieved  in  [3]  through  a  rather 
elaborate  sequence  of  operations  on  the  Taylor  series  coefficients  of  A(e), 
involving  cascaded  projections  onto  progressively  slower  subspaces.  A  major 
advantage  of  this  result  is  its  generality  —  with  it  we  can  analyze  general 
systems  as  in  (1.1)  without  assuming  some  special  form.  A  price  that  is  paid 
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for  this,  however,  is  that  the  results  and  procedures  developed  are  rather 
complicated,  involve  the  computation  of  numerous  pseudo-inverses,  and 
generally  do  not  lend  themselves  to  easy  interpretation  or  computation. 

The  work  presented  in  this  paper  bridges  the  gap  between  the  intuitively 
and  computationally  simple  but  somewhat  restrictive  results  of  [1],  [2]  and 
the  quite  general  but  rather  complicated  ones  in  [3],  The  key  to  constructing 
this  bridge  is  an  examination  of  the  algebraic  structure  of  A(e)  considered  as 
a  matrix  over  the  ring  W  of  functions  of  e  that  are  analytic  at  e  =  0.  In 
particular,  by  considering  the  Smith  form  of  A(e)  we  not  only  provide  the 
basis  for  transforming  a  general  system  (1.1)  to  its  explicit  form,  but  also 
make  clear  the  role  of  the  invariant  factors  of  A(e)  in  specifying  the  time 
scales  present  in  the  dynamics  (1.1),  a  role  that  is  suggested  but  not 
developed  in  [3].  This  approach  provides  some  valuable  additional 
perspectives  on  the  results  in  [1]  -  [3],  and  it  also  allows  us  to  consider 
and  solve  a  number  of  additional  problems.  Several  of  these  are  presented  in 
the  later  sections  of  this  paper,  while  others  will  be  the  subject  of  future 
papers.  We  note  here  that  another  approach  to  the  main  results  of  [3]  is 
described  in  [6],  which  proceeds  by  transforming  A(e)  to  a  block-diagonal  form 
that  is  similar  to  it.  There  is  a  clear  point  of  contact  between  our  work  and 
the  results  in  [6],  as  our  proof  in  Section  4  of  the  sufficiency  of  certain 
conditions  for  the  existence  of  a  time  scale  approximation  is  much  in  the 
spirit  of  the  methods  in  [1],  [2],  and  [6].  On  the  other  hand,  our  results  go 
significantly  farther  than  previous  efforts  in  that  we,  for  the  first  time, 
make  clear  the  role  of  the  Smith  form  and  the  invariant  factors  of  A(e)  and 
present  a  procedure  that  minimizes  the  number  of  e-dependent  computations 


required  by  identifying  and  discarding  non-critical  t-dependent  terms  in  A(e) 
and  in  its  explicit  form. 

In  the  next  section  we  introduce  a  new  definition  of  what  we  call  a 
strong  time  scale  decomposition.  Based  on  this,  we  present  a  new  result  that 
allows  us  to  state  a  strengthened  version  of  the  main  result  in  [3]  and  to 
obtain  a  criterion  for  identifying  higher-order  terms  in  a  system  matrix  A(e) 
that  can  be  discarded  without  affecting  the  investigation  of  the  existence  of 
strong  time  scale  behavior.  In  Section  3  we  then  introduce  the  Smith  form  of 
A(e)  and  use  it  to  transform  (1.1)  to  explicit  form.  We  also  perform  some 
initial  analysis  that  allows  us  to  focus  subsequent  discussions  on  the  case  in 
which  A(e)  is  Hurwitz  for  0  <  e  <  eQ  for  some  >  0.  In  Section  4  we  develop 
what  can  be  viewed  as  a  generalization  of  the  procedure  in  [1],  [2]  to  analyze 
systems  in  explicit  form.  This  produces  both  a  set  of  necessary  and 
sufficient  conditions  for  a  system  to  have  a  strong  time  scale  decomposition 
and  a  procedure  for  constructing  the  corresponding  strong  multiple  time  scale 
approximation. 

With  these  results  established,  we  can  then  consider  two  important 
extensions.  In  Section  5  we  consider  a  generalization  of  the  definition  of  a 
time  scale  approximation  that  allows  us  to  construct  such  approximations  for  a 
large  class  of  systems  violating  the  conditions  of  Section  4.  In  Section  6  we 
address  the  problem  of  modifying  and  controlling  the  time  scales  of  the  system 
*(t)  =  A(e)x( t )  +  B(e)u( t)  (1.11) 

through  the  use  of  feedback 


(|  *•* A  A  kV|^' La’Wj.^Lt'tVLi ‘ t >t! f ■  LWiU il.Wil»W il.T* 
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2.  Well-Defined  Multiple  Tiae  Scale  Behavior 

To  begin  this  section  we  give  two  different  definitions  of  what  one  might 
mean  by  well-defined  multiple  time  scale  behavior.  The  first  of  these  is 
essentially  the  standard  definition  that  is  stated  or  implied  in  previous 
treatements.  The  seond,  stronger  definition  is  new,  as  it  requires  the 
consideration  of  an  entire  family  of  systems.  By  introducing  this  definition 
we  can  make  several  new  observations  concerning  time  scale  decompositions  and 
can  give  a  stronger  interpretation  of  the  results  in  [3]. 


Def ini tion  2.1*  The  system  (1.1)  has  a  multiple  time  scale 

decomposition  if  there  exist  constant  matrices  At ,  A~ .  A  ,  T 

12  n 

and  integers  0  <  k,  <  k0  <...<  k  such  that 
-  l  2  n 

k  k  k 

lim  sup  ||eA^e^t  -  T  *exp  (diag  e  2A„ . e  nA  t  }T  |  |  =  0 

elO  t>0  1  ^  nJ  % 

(2.1) 

In  this  case  we  say  that  [{A^ ,  {k^  ,T]  defines  a  multiple  time 
scale  decomposition  of  (1.1)  or  of  A(e) . 


To  introduce  the  second  definition  we  first  need  the  following:' 


Definition  2.2:  The  perturbed  family  ?{A(e)}  associated  with  the 
matrix  A(e)  is  defined  as  follows: 

*{A(e)}  =  <U(e)A(e)V(e)|U(0)  =  V(0)  =  I)  (2.2) 

Definition  2.3:  The  system  (1.1)  has  a  strong  multiple  time  scale 

decomposition  if  there  exist  constant  matrices  A. ,  A0, ....  A  ,  T 

12  n 


Throughout  this  paper  we  assume  that  all  matrix  functions  of  t  are  analytic 
at  zero. 
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and  integers  0  <  <...<knsuch  that 


lim  sup  | |e^^£^ t-T  *exp  {diag 
eiO  t>0 


k  k 

e  *Aj . e  nA 


for  al  1  F(e)  €  ?{A(e.)}.  In  this  case  we  say 
defines  a  strong  time  scale  decomposition  of 


n  t }T | |  =  0  (2.3) 

that  [{A.}.{k.},T] 
(1.1)  or  of  A(e) . 


Clearly  the  second  of  these  definitions  is  significantly  stronger  that 

the  first.  Intuitively  the  elements  of  5{A(e)}  should  be  thought  of  as  mild 

perturbations  of  A(e),  and  the  strong-sense  definition  requires  that  any  such 

perturbation  must  result  in  a  system  that  has  the  same  time  scale 

decomposition  as  (1.1).  More  precisely,  an  immediate  consequence  of 

Definition  2.3  is  that  if  A(e)  has  a  strong  time  scale  decomposition,  then  any 

G(e)  €  ?{A(e)}  is  asymptotical ly  equivalent  to  A(e),  i.e. 

lim  sup  ||eA^£^t  -  eG^t||  =  0  (2.4) 

eiO  t>0 

To  illustrate  these  ideas  let  us  consider  several  examples.  First,  note 
that  the  scalar  system 

x(t)  =  x(t)  (2.5) 

trivially  has  a  time  scale  decomposition  according  to  Definition  2.1.  but  not 
according  to  Definition  2.3  since  (1+e)  €  ?{1)  is  not  asymptotically 
equivalent  to  1 .  On  the  other  hand,  it  is  not  difficult  to  check  (and  is 
immediate  from  the  results  in  several  papers)  that 

x(t)  =  -  x(t)  (2.6) 

does  have  a  strong  time  scale  decomposition. 

Consider  next  the  system  matrix 
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o  r 

A(e)  =  (2.7) 

.-i  0. 

This  matrix  has  a  trivial  time  scale  decomposition,  but  it  does  not  have  a 
strong  time  scale  decomposition,  since  it  is  not  asymptotically  equivalent  to 
the  matrix 


(2.8) 


Finally,  we  note  that 


0  f 

A(e)  =  (2.9) 

lo  oj 


does  not  have  a  strong  time  scale  decomposition  since  it  is  not  asymptotically 
equivalent  to 


These  examples  indicate  that  there  are  problems  when  there  are 
eigenvalues  that  are  in  the  right-half  plane,  are  purely  imaginary,  or  are 
zero  with  nontrivial  Jordan  blocks.  To  see  that  these  examples  span  all 
possible  cases,  we  need  to  briefly  re-examine  and  strengthen  the  main  result 
in  [3],  In  particular,  although  it  is  not  discussed,  [3]  in  fact  provides  the 
basis  for  determining  if  a  system  has  a  strong  time  scale  decomposition  and 
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for  constructing  that  decomposition.  The  system  considered  in  [3]  is  the 
singularly  perturbed  LTI  system  (1.1)  with  a  slight  change  in  notation  whose 
purpose  will  become  clear  shortly. 

*(t)  =  AQ(e)x(t)  (2.11) 

where  the  matrix  A^(e)  is  an  analytic  function  of  t  at  t  =  0.  Suppose  A^(e) 

has  eigenvalues  A^e) . A^(e)  where  A^(t)  -*  0.  t  -*  0,  i  =  1 . m  <  n. 

Then  the  total  projection  for  the  zero-group  of  eigenvalues  of  A^(e).  P^ft). 
is  the  projectior;  onto  the  subspace  spanned  by  eigenvectors  and  generalized 


eigenvectors  corresponding  to  A.(t) . A^(t)  of  A^(e)  [^]- 

Since  A„(fc)  is  analytic  at  e=0,  it  has  a  series  expansion  of  the  for 


p=0 


(2.12) 


It  can  be  proven  [3.  7]  that  if  has  semisimple  nu 1  I s t rue ture  (SSNS)  — 
i.e.  if  its  zero  eigenvalue  is  semisimple,  that  is.  lias  geometric  multiplicity 
equal  to  its  algebraic  multiplicity  —  then  the  matrix 


Al(<0  =  P0(e)A0(e)/e  (2.13) 

has  a  series  expansion  of  the  form 

00 

Aj(e)  =  2  ePF  (2.14) 

p=0  y 

(otherwise  A^(e)  will  have  e  ^  terms).  If  F^q  also  has  SSNS  we  define  Agfc) 
as 

*2(t)  =  PjttjAjfeJ/t  =  PjleJPQfeJA^el/e2 

■  jo  **  (215) 

where  Pj(e)  is  the  total  projection  for  the  zero-group  of  eigenvalues  of 
A^e).  This  process  can  be  continued  until  it  terminates  at 


if  the  matrix  F  does  not  have  SSNS  or  if  rankF„ _+rankF +rankF  _  equals 
riO  00  10  nO 


i  he 


normal  rank  of  A()(f  ).  i.e.  t Fie  constant  rank  that  Aq(&)  takes  on  some 
interval  (0.  t(J],  A  matrix  AQ(t)  is  said  to  satisfy  the  multiple  semisimple 
nu 11  st  rue  t  ure  (MSSNS)  condition  if  the  latter  of  these  conditions  holds.  If 


in  .addition,  all  F^q  are  semi  stable 


i.e.  if  for  each  k  F,  _  has  SSNS  and  all 

kO 


of  its  nonzero  eigenvalues  have  strictly  negative  real  pans  —  then  we  say 
that  Aq(c)  satisfies  the  multiple  sem i s tab i 1 i t y  (MSST)  condition. 

The  main  result  of  [.')]  is  that  if  A(t)  satisfies  MSST,  then  (i) 


kO 


T  diag(0 . 0.  A.  ,0. 


,  0)T 

otherwi se 


.  k  =  k. 


(2.17) 


for  some  nonsingular  T,  semistable  A^  ,  and  uniquely  determined  integers  k.; 
and  (ii)  A(t)  has  a  time  scale  decomposition  in  the  sense  of  Definition  2.1. 

On  the  other  hand,  as  our  examples  (2.3).  (2.7),  (2.9)  show.  MSST  is  not 
necessary  for  A(e)  to  have  a  time  scale  decomposition.  What  we  show  in 
Theorem  2.1  is  that  MSST  is  necessary  and  sufficient  for  A(e)  to  have  a  strong 
time  scale  decomposition. 

In  order  to  prove  our  strengthened  version  of  the  main  result  in  [3]  we 
need  two  results. 


Proposi t ion  2.1-  Let  G(t)  f  ?{A(c)).  Then 


Vk 


(2. IS) 


where  the  superscripts  "G"  and  "A"  denote  the  sequences  defined  in 
(2.12)  -  (2.16)  for  G(c)  and  A(t),  respectively 
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Proof •  See  Appendix  A. 

Proposition  2.2:  Suppose  that  [{A^ } , (k^ } ,T]  defines  a  multiple 
time  scale  decomposition  of  Aq(&)  and  suppose  further  that 

Aj . A^  are  semistable.  Then  (2.17)  holds  and  hence  A(e) 

satisfies  the  MSST  condition. 

Proof :  See  Appendix  B 

We  can  now  state: 

Theorem  2.1:  The  system  (2.11)  has  a  strong  time  scale 
decomposition  if  and  only  if  Aq(&)  satisfies  the  MSST  condition. 

Proof :  As  stated  previously,  it  is  proved  in  [3]  that  the  MSST  condition  is 

sufficient  to  satisfy  the  weaker  Definition  2.1.  That  this  condition  is 

sufficient  for  the  stronger  definition  follows  directly  from  Proposition  2.1 

The  proof  of  necessity  is  also  straightforward.  Specifically  if  A^fe)  has  a 

strong  tune  scale  decomposition  as  in  (2.1),  then,  thanks  to  Proposition  2.2 

all  we  need  to  show  is  that  the  A^  must  be  semistable.  This  can  be  done  by 

contradiction.  Specifically,  if  Ai  is  not  semistable,  then  it  has  a 

right-half  plane  eigenvalue,  a  pair  of  purely  imaginary  eigenvalues,  or  a 

nontrivial  Jordan  block  corresponding  to  the  0  eigenvalue.  Showing  that  any 

of  these  conditions  preclude  the  existence  of  a  strong  time  scale 

decomposition  is  a  minor  variation  on  our  previous  discussion  of  the  three 

examples  (2.5),  (2.7),  (2.9). 

Note  that  if  A(e)  is  invertible  for  e  €  (0,  £q],  the  A^  in  the  strong 
time  scale  decomposition  are  all  Hurwitz. 


Finally,  it  is  also  shown  in  [3]  that  if  AQ(e)  does  not  satisfy  MSST, 
then  for  some  q  the  limit  as  eiO  of 

exp  [AQ(e)t/eC1]  (2.19) 

does  not  exist.  This  indicates  that  a  failure  of  the  MSST  condition  does 
correspond  to  some  type  of  noncovergent  behavior.  However,  the  precise 
meaning  and  interpretation  of  this  could  not  easily  be  exposed  without  the 
concept  of  a  strong  time  scale  decomposition.  Indeed,  in  addition  to 
providing  us  with  Theorem  2.1,  this  machinery  makes  it  far  simpler  to  prove 
the  nonexistence  of  the  limit  of  (2.19).  Furthermore,  we  now  see  that  to 
verify  the  MSST  condition  and  to  construct  a  time  scale  decomposition  for 
A(e),  we  can  equivalently  examine  these  questions  using  any  element  of  ?{A(e)} 
—  i.e.  any  such  element  must  generate  the  same  sequence  F^q  if  a  strong  time 
scale  decomposition  exists.  Of  course  we  can  equivalently  consider  any 
element  of  ?{SA(e)S  *}  where  S  is  any  constant  invertible  matrix.  We  make 
use  of  these  facts  in  the  next  section  to  transform  an  arbitrary  A(e)  to  its 


explicit  form. 
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3.  Explicit  Form 

As  mentioned  in  Section  1,  our  new  approach  employs  the  Smith 
decomposition  of  A(e)  over  the  ring  W  of  functions  of  £  that  are  analytic 
at  e  =  0  (see  [4],  [5]).  The  units  of  W  are  elements  that  do  not  vanish  at 
e  =  0.  That  is,  since  any  element  of  W  can  be  expanded  in  a  Taylor  series. 


a0  +  ai£  +  a2e  +  • • • 


(3.1) 


we  see  that  the  set  of  units  are  those  elements  with  a^  ^  0.  It  is  also 
easily  seen  that  W  is  a  Euclidean  domain,  with  the  degree  or  order,  0(d(e)), 
of  any  element  d(e)  €  W  being  defined  as  the  order  of  the  first  nonzero  term 
in  its  Taylor  expansion.  Therefore  A(e)  has  a  Smith  decomposition 

A(£)  =  P(£)D(£)Q(£)  (3.2) 

where  P(e)  and  Q(fc)  are  unimodular,  i.e.  jP(e)  |  and  |Q(e) j  are  units  (and  thus 
P  *(e)  and  Q  * (e )  are  matrices  over  W)  or,  equivalently 


|P(0) |  *  0  .  |Q(0) |  *  0 


(3.3) 


and 

k  k 

D( £ )  =  diag  ( £  lI . e  nI .  0)  (3.4) 

where  0  <  kj  <  <  ...  <  kn  are  integers,  the  identity  matrices  I  may  have 

different  dimensions,  and  the  0  matrix  is  only  present  if  A(e)  is  singular  in 

k . 

a  neighborhood  of  £  =  0.  The  e  are  called  the  invariant  factors  of  A(e) 
Actual  computation  of  such  Smith  decompositions  is  discussed  in  [4]  and  [5] 

(in  the  terminology  of  [5],  what  is  required  is  to  transform  A(e)  to  the 
matrix  D(e)Q(e)  which  is  "row-reduced  at  0"  through  row  operations  embodied  in 
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P  1(e)).  Without  loss  ul  generality  we  assume  from  here  on  that  =  0;  this 
can  always  be  obtained  by  a  change  of  time  scale  in  (1.1). 

Rather  than  working  with  the  system  (1.1),  we  consider  an  e-independent 
change  of  variables 

y(t)  =  P_1(0)x(t)  (3.5) 


so  that 

y(t)  =  P_1 (0)P(e)D(e)Q(e)P(0)  (3.6) 

Next  we  note  that  if  we  define  the  constant  matrix 

A  =  Q(0)P(0)  (3.7) 

then 

D(e)A  €  ?{P_1 (0)P(e)D(e)Q(e)P(0) }  (3.8) 

(premultiply  P  * (0)P(e)D(e)Q(e)P(0)  by  P  ^(e)P(O)  and  postmultiply  by 
P  1(0)Q  ^(e)Q(0)P(0) ) .  Therefore,  we  arrive  at  the  explicit  form  of  (1.1): 

z  =  D(e)Az  (3.9) 

which,  if  we  express  A  in  block  form  with  blocks  compatible  with  those  in 
(3.4),  can  be  written  as 


’*1 

A11 

A12 

Al,n+1 

Z1 

= 

k0 

k„ 

k„ 

*2 

2  A 

.  A2i 

,  2  . 

.  A22 

2 

6  A2 , n+ 1 

Z2 

• 

k 

k 

k 

• 

z 

n 

n  A 

6  A  , 
nl 

11  A 

6  An2  • 

n  A 

e  A  . . 

n,n+l 

z 

n 

.  Vi. 

0 

0 

0 

.  Zn+1. 

(3.10) 


j 


Let  us  make  several  comments  about  the  previous  transformations.  First 


of  all,  note  that  every  element  of  ?  (A(e)}  has  the  same  explicit  form. 
Secondly,  if  A(e)  does  not  have  a  strong  time  scale  decomposition,  then,  as  in 
the  examples  in  Section  2,  there  is  no  reason  to  expect  that  (3.9)  is  a  good 
approximation  of  (3.6)  (and  therefore  of  (1.1))  in  that  the  two  systems  need 
not  be  asymptotically  equivalent.  However,  if  any  of  the  systems  (1.1), 

(3.6),  or  (3.9)  has  a  strong  time  scale  decomposition,  then  they  all  do,  and 
(3.9)  is  asymptotically  equivalent  to  (3.6).  Therefore,  we  can  focus  on  the 
explicit  form  if  we  are  interested  in  strong  time  scale  decompositions. 
Finally,  note  that  the  system  (3.10)  is  an  obvious  generalization  of  (1.2), 
and  this  observation  provides  the  basis  for  our  development  in  the  next 
section.  Before  doing  this,  however,  we  first  conclude  this  section  by 
showing  how  we  can  deal  with  the  0  diagonal  block  in  D(e)  so  that  hereafter  we 
can  focus  attention  on  the  case  in  which  there  is  no  such  block,  i.e.  the  case 
in  which  A(e)  is  Hurwitz  for  e  €  (0,  e^] . 

Specifically,  let  us  write  D(e)  in  (3.4)  as 

D(fc)  =  diag  (D^e).  0)  (3.11) 

(so  that  Dj(e)  consists  of  all  of  the  nonzero  invariant  factors),  and  let  us 
express  A  in  (3.7)  in  blocks  compatible  with  (3.11). 

G11  G12 

A  = 

G21  G22 
We  then  have  that 

D(e)A  =  (3.13) 

10  0 

Note  that  (G^  G^)  *1as  row  ran^  since  A  is  invertible.  In  fact,  it  is 

immediate  from  the  development  in  the  next  section  that  D(e)A  has  MSSNS  only 


(3.12) 
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1£  Gji  is  invertible.  Therefore,  as  a  first  step  in  our  overall  procedure,  we 
check  the  inver tibi  1  i ty  of  G^.  If  it  is  not  invertible,  then  we  immediately 
know  that  (3.9)  and  hence  (1.1)  do  not  have  strong  time  scale  decompositions. 
If  Gjj  is  invertible,  we  perform  the  following  e-independent  transformation  of 


(3.9) 


G11  1C121 


(3.14) 


so  that 


fDl<‘>GU  °1 


(3.15) 


From  this  point  on  we  can  focus  completely  on  the  lower-dimensional,  explicit 
form  matrix  D^(e)G^  which  is  invertible  for  e  €  (0,  £q] •  If  this  has  a 
strong  time  scale  decomposition,  then  so  do  (3.9)  and  (1.1),  and  the 
construction  of  the  time  scale  approximations  for  these  systems  from  the  one 
for  (3.15)  involves  the  obvious  reversal  of  the  steps  taken  to  obtain  (3.15) 
f rom  (1.1). 


4.  Strong  Multiple  Tire  Scale  Decompositions  of  Systems  in  Explicit  Form 


Based  on  the  development  and  discussion  in  the  previous  section,  we  now 
focus  attention  on  the  following  system  in  explicit  form 

z(t)  =  D(e)Az(t)  (4.1) 

k2  kn 

where  D(e)  =  diag  (I,  e  I . £  I)  and 


A11 

A12  ‘  ‘ 

■  A,  1 
In 

*n:j! 

A 

A 

A„ 

_ 

*21 

22 

2n 

w1  :  z 

v 

<3 

■  < 

.  A 

ruij 

(4.2) 


is  invertible.  The  reasons  for  the  notation  introduced  in  (4.2)  will  become 
clear  shortly  (here  the  dashed  line  in  both  matrices  are  in  the  same 
locations,  so  that  =  A^,  =  [A^-'-A^]-  etc)’ 

One  direct  approach  to  determining  necessary  and  sufficient  conditions 
under  which  (4.1)  (and  thus  (1.1))  has  a  strong  time  scale  decomposition  is  to 
identify  explicitly  the  projections  and  similarity  transformations  used  in  [3] 
to  check  for  MSST  and  to  ■  Stain  the  multiple  time  scale  decomposition 
described  in  Theorem  2.1.  This  is  done  in  detail  in  [8].  What  we  do  in  this 
section  is  to  follow  an  approach  that  makes  use  of  the  results  in  Section  2  to 
obtain  a  set  of  necessary  and  sufficient  conditions  and  a  procedure  for 
constructing  a  multiple  time  scale  decomposition  that  is  much  more  in  the 
spirit  of  [1]  and  [2].  Based  on  our  initial  review  of  the  analysis  of  (1.2), 
it  should  not  come  as  a  surprise  that  successive  Schur  complements  of  A  play 
am  important  role  in  our  development.  Also,  since  we  are  focusing  on  strong 
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time  scale  decompositions,  we  have  the  luxury  of  throwing  away  many  of  the 
e-dependent  terms  that  arise  as  we  proceed.  Specifically,  whenever  we  run 
into  a  unimodular  matrix  U(e)  multiplying  our  system  matrix  on  the  left  or 
right,  we  can  replace  i t  by  U(0)  and  continue.  Either  both  of  these  systems 
have  the  same  strong  time  scale  decompositions  or  nei ther  one  has  such  a 
decomposi tion. 

The  basic  idea  behind  the  approach  we  use  here  is  to  block-diagonal ize 
D(e)A.  We  do  this  in  stages,  "peeling  off"  one  time  scale  of  (4.1)  at  a  time, 
starting  with  the  fastest.  To  begin,  let  us  introduce  some  notation. 
Specifically,  let  D^(e)  =  D(e),  =  A,  and 


D(e)A  =  D1(e)A1  = 


k2 

U  2S,(e) 


(4.3) 


Sx{e)  =  D2(e)W1 
F^e)  =  D2(e)Z1 

k-Tk9  krTk9 

D2(e)  =  diag(  I ,  e  13  Z1 . £  n  ZI) 


(4.4) 

(4.5) 


(4.6) 


(here  the  dimensions  of  the  (n-1)  identity  matrices  in  (4.6)  are  the  same  as 

the  last  (n-1)  blocks  in  D(e)). 

As  a  next  step  we  prove  the  following: 

Lemma  4.1:  Consider  the  constant  matrix 
fN  Ll 


0  0 


(4.7) 


where  (N,  L)  has  full  row  rank  and  N  is  square.  Then  M  has  SSNS  if  and  only 
if  N  is  invertible. 
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Proof :  Suppose  N  is  invertible.  Then 


-1  1 
I  N  L 

M 

I 

►J 

i 

z; 

i 

N 

o' 

0  I 

0 

I  J 

0 

0. 

M  =  (4-8) 

[0  I  J  [0  I  J  [0  0  J 

which  clearly  has  SSNS.  On  the  other  hand,  if  N  is  not  invertible,  then  there 
exists  x  ^  0  so  that  Nx  =  0.  Furthermore,  since  (N,  L)  has  full  row  rank,  we 
can  find  Xj  and  so  that  Nx1  +  Lxg  =  x.  If  we  then  define 

hi 

z  =  1  (4.9  ) 

L  2J 

we  have  that  Mz  /  0  but  M^z  =  0,  showing  that  M  does  not  have  SSNS. 


Letting  t  =  0  in  (4.3),  we  have 


D1(0)A1  = 


hi  R1 


0  0 


(4.10) 


Since  is  invertible.  [A^,  Rj]  has  full  row  rank.  Consequently,  from  Lemma 
4.1  we  see  that  the  system  matrix  (4.10)  describing  evolution  at  the  fastest 

IV  /\> 

time  scale  has  SSNS  if  and  only  if  A^  is  invertible.  Suppose,  then,  that  A^ 
is  invertible.  Consider  the  similarity  transformation 


G(e)  = 


'I 

~  -i  - 
A11  R1 

D1(e)A1 

I 

-*nRi 

P 

I 

0 

I 

Ah  +  6  ^iflRiD2(*)Wi 


e\{e)k2 


(4.11) 


where 


*2  *  Z1  -  Vll  R1 


(4.12) 


which  is  invertible  (since  and  in  (4.2)  are  both  invertible).  Note 


further  that 


G(fc)  =  U(e) 


2D2(e)A2J 


(4.13) 


where  U(0)  =  V(0)  =  I  (see  Appendix  C) . 

Since  we  are  interested  in  strong  time  scale  decompositions,  we  can 
discard  U(e)  and  V(e).  From  Proposition  2.2  and  Theorem  2.1  we  can 

^  _ 

immed iate ly  conclude  that  for  diag  (A^ ,  a  u^{a)k^)  to  have  a  strong  time 
scale  decomposition,  A^  must  be  Hurwitz.  Furthermore,  we  have  now  reduced 
the  problem  to  the  examination  of  the  explicit  form  matrix  Dg^JAg  with  one 
fewer  time  scale. 

Consider  now  the  following  recursion  beginning  with  A^  in  (4.2)  and 
defined  recursively  as  follows 


[A.  .  R. 

n  l 


Ai+1  =  zi  -  Vii  Ri 


(4.14) 


(4.15) 


Here  the  block  size  of  each  A^  is  the  same  as  that  of  the  ith  block  in  the 

g 

original  explicit  form  systems  (4.1),  (4.2)  .  Using  the  results  of  Section  2 


Note  that  at  the  last  step  A  =  A 

n  nn 


then  yields  the  following: 


Theorem  4.1:  The  explicit  form  system  (4.2)  has  MSSNS  if  and  only 

if  each  A. .  is  invertible.  Furthermore,  the  system  (4.2)  satisfies 
n 

the  MSST  condition,  and  hence  has  a  strong  time  scale. decomposi t ion 
if  and  only  if  each  of  the  A^  is  Hurwitz.  In  this  case 


D(e)At  _1  ^  k~,  k 

lim  sup  |  |e  -  T  exp  {diagEA^.e  . e  A^jtJTH  =  0 

eiO  t>0 


where 


IP  I  J 


I 

0 

Ti  = 

0 

~  -1  -1 
1  Aii  Ri 

0  I 

(4.16) 

(4.17) 

(4.18) 

i  >  1  (4.19) 


(here  the  upper  left-hand  identity  block  is  of  dimension  equal  to 
the  first  ( i— 1 )  blocks  of  (4.1),  (4.2)). 


We  close  this  section  with  several  final  comments.  First,  note  that  the 
recursive  procedure  just  described  for  peeling  off  successively  slower  time 
scales  actually  yields  a  sequence  of  approximations  over  successively  longer 
time  intervals,  i.e. 

, .  II  D(e)At  ,  ^  k  ^ 

o  1,6  -T  ^xp  (diag[An . e  rArr.O....O]t)T||  =  0 

te[0,e  r+1) 


(4.20) 


(see  [3,  6]  for  similar  comments).  Secondly,  we  note  that  an  alternate 

approach  to  showing  the  sufficiency  of  the  conditions  in  Theorem  4.1  is 

presented  in  [9]  using  an  approach  much  in  the  spirit  of  [6].  Specifically. 

consider  the  following  equations 
^  k 

Rj  +  AjjL^e)  -  e  ^(e^F^e)  +  S^eJL^e)]  =  0  (4.21 


H1(e)[An  -  e  ^(eJS^e)]  -  e  ^(e)  +  +  S^e)  =  0 

(4.22 

It  is  straightforward  to  check  that  these  equations  have  solutions  L^(e)  and 
H^(e)  for  e  small  enough,  that 

1^(0)  =  -  ^n~\.  ^(0)  =  -  S^O)^"1  (4.23 

that 


I 


^(e)  = 


~L1(fe) 


I-e  ^(eJLjU) 


(4.24 


is  unimodular,  and  that  the  similarity  transformation  specified  by  T^(e) 
block-diagonal izes  D^(e)A^,  i.e. 


T1(£)D1(e)^1T~1(fc)  = 


'Gj(e)  0 

0  6  G2(6)J 


(4.25 


where 
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Gj (e)  =  An-e  ^(eJS^e)  G2(d)=F1(e)+S1  (04(e)  (4.26) 

/v  ___ 

Noting  then  that  G^(O)  =  and  that  G^Cfc)  =  D2(e)A2C(e)  where  C(0)  =  I.  we 
can  conclude  that  D^feJAj  has  a  strong  time  scale  decomposition  if  and  only  if 


T i (0)  (4.27) 


does,  where 


T1  \0)  = 


lll 


0 


0 


(e)Ar 


Tj  *(0) 


'I 

0 


I 


(4.28) 


This  process  can  then  be  iterated  to  consider  the  next  time  scale. 

Comparing  this  procedure  to  that  described  previously,  and  in  particular 
to  (4.11)  and  the  subsequent  development,  we  see  that,  thanks  to  Theorem  2.1, 
we  do  not  have  to  do  quite  so  much  work  (although,  as  described  in  Appendix  A, 
we  actually  use  this  full  block-diagonal ization  procedure  in  the  proof  of 
Proposition  2.1).  Rather,  instead  of  fully  block-diagonalizing  D^(e)A^  using 
the  full  Tj(e),  we  simply  use  T^(0),  the  key  being  that  we  have  raised  the 
order  of  the  upper  right-hand  element  of  (4.11)  suf f iciently  so  that  (4.13) 
holds.  In  a  sense  what  we  have  done  in  (4.11)  is  a  first  step  in  an  itein  ivc 
approach  to  block-diagonalizing  4(6)4.  Specifically,  think  of  the 
transformation  in  (4.11)  as  an  attempt  to  approximately  null  out  the  (1,  2) 
block  of  D^eJAj  by  raising  its  order.  If  we  then  attempt  to  approximately 


null  out  the  (2,  1)  block  of  G(e)  (using  a  lower-block- 1 r iangu lar  similar  if  , 
transformation),  we  will  raise  the  order  of  this  term.  Carrying  this  process 
on  for  a  number  of  steps  we  obtain  better  and  better  approx i mn i e  block 
diagonalizations  and  hence  have  a  series  expansion  for  T  (<  )  What  we  have 
shown  here  is  that  when  looking  for  strong  time  scale  decompositions,  we  can 
stop  after  the  first  term  in  the  series.  in  the  next  section  we  describe  a 
procedure  for  constructing  a  weaker  form  of  a  time  scale  decomposition  for 
systems  not  satisfying  the  MSST  condition.  This  procedure  requires  keeping 
additional  terms  of  the  series  or,  equivalently,  performing  the  iterative, 
approximate  block-diagonal izat ion  procedure  for  more  than  one  iteration. 


5.  Tine  Scale  Decompositions  for  Systems  Without  MSST 


In  this  section  we  describe  a  procedure  for  constructing  a  somewhat 

weaker  time  scale  decomposition  for  systems  that  do  not  satisfy  the  MSST 

condition.  To  motivate  and  illustrate  the  essential  ideas  behind  this 

procedure,  we  begin  with  an  example.  Specifically,  consider  the  system  matrix 
f-e  ll 


AO)  = 


(5.1) 


l-l  -£J 

Since  A(0)  is  not  semistable  we  immediately  see  that  this  matrix  does  not  have 
a  strong  time  scale  decomposition.  In  fact,  ■:  t  is  not  difficult  to  see  that 
it  does  not  even  have  a  time  scale  decomposition  in  the  sense  of  Definition 
2.1.  The  reason  for  this  stems  from  the  requirement  that  the  system  matrices 
Aj,  .  in  (2.1)  be  independent  of  e.  Examining  A(e)  in  (5.1)  we  see  that 

its  eigenvalues  (-  e  ±  j)  have  the  property  that  their  real  parts  are  of 
higher  order  in  e  than  their  imaginary  parts.  Consequently,  when  we  attempt 
to  use  a  constant  system  matrix  to  approximate  (5.1)  we  throw  away  the  crucial 
damping.  From  this  perspective  it  seems  evident  that  what  one  should  seek  to 
do  in  this  case  is  to  keep  at  least  some  of  the  e-dependent  terms  in  A(e)  in 
order  to  preserve  its  principal  damping  characteristics.  The  procedure  we 
develop  in  this  section  does  exactly  that. 

We  begin  our  development  with  the  following 


Def ini tion  5.1:  Let  A(e)  be  Hurwitz  for  e  €  (0,  and  let  the 
Smith  form  of  A(e)  be  as  in  (3.2)  with  D(e)  = 
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lim  sup|  |eA^fc^t-T  ^xpfdiagO-te ) . tA  (e) . en  1 A  ( ) ] t }T  |  |  =  0 

eiO  t>0  u  l  n 


(5.2) 


where  T  is  a  constant  matrix  and  each  of  the  A.(e)  has  the 

properties  that  A^O)  is  invertible  and  each  of  its  purely 

imaginary  eigenvalues  is  semisimple  (i.e.  has  algebraic 
multiplicity  equal  to  its  geometric  multiplicity). 


Let  us  make  several  comments  about  this  definition.  First,  using  the 
procedure  described  at  the  end  of  Section  3  we  can  actually  weaken  the 
assumption  that  A(e)  is  Hurwitz  by  assuming  only  that  A(e)  is  semistable  for 
e.  €  (0,  £q]  (so  that  there  may  be  a  0  block  in  D(t));  however  for  simplicity 
here  we  use  the  stronger  assumption.  Also,  the  assumption  that  D(e)  has  the 
particular  form  stated  in  the  definition  is  no  real  restriction  and  again  we 
include  it  here  for  convenience  only  (if  some  power  of  £  between  0  and  n-1  is 
not  an  invariant  factor,  then  the  corresponding  step  of  our  procedure  is 
simply  dropped).  Finally,  let  us  discuss  the  assumptions  on  A^O).  Note 
first  that  requiring  A^fO)  to  be  invertible  is  equivalent  to  assuming  that 
A(e)  has  MSSNS,  while  the  further  semisimplicity  assumption  eliminates 
matrices  such  as 


-£ 

1 

1 

0 

-1 

-£ 

0 

1 

0 

0 

-£ 

1 

0 

0 

-1 

-£ 

A  A  A-V  AAA  A  -■  A  -VA-S  -VIN  A  -VA-’, 


which  are  Hurwitz  for  e  >  0  but  for  which 

sup  | |exp{F(e) t}  |  | 
t>° 

grows  without  bounds  as  e  i  0.  In  essence  what  we  are  considering  in  this 
section  is  the  extension  of  our  theory  of  time  scale  decompositions  to  include 
A(e)‘s  with  eigenvalues  that  converge  to  points  on  the  imaginary  axis  other 
than  the  origin.  Consequently,  it  is  not  surprising  that  the  multiple 
semisimplicity  condition  is  extended  to  include  all  eigenvalues  converging  to 
the  imaginary  axis. 

Definition  5.2 :  A  matrix  A(e)  has  multiple  semisimple  imaginary 
e i gens t rue ture  (MSSIES)  if  it  has  MSSNS  and  if  each  of  the  purely 

imaginary  eigenvalues  of  each  of  the  A_  defined  in  Theorem  4.1  is 

semi simple . 

Essentially  by  definition  we  have  that  MSSIES  is  necessary  for  A(e)  to  have  a 
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weak  time  scale  decomposition.  In  fact,  the  procedure  we  describe  in  this 
section  proves  the  following: 

Theorem  5.1:  Let  A(e)  be  Hurwitz  for  e  €  (0.  e^].  Then  A(e)  has  a 

weak  multiple  time  scale  decomposition  if  and  only  if  it  has 
MSSIES. 


Indeed  if  this  is  not  the  case,  then  (5.2)  leads  to  a  contradiction,  since 

lim  sup  ||e^£^t||  =  ® 
tlO  t>0 

but  exp{diag[A(-)(£) . en  ^A^ _^(e)]t}  is  uniformly  bounded. 


TTT  YTT 


For  the  remainder  of  this  section  we  assume  that  A(e)  is  Hurwitz  and  has 
MSSIES.  As  a  first  step  in  our  procedure,  we  transform  the  dynamics  of  (1.1) 
in  a  manner  similar  to  that  used  in  Section  2.  Specifically,  let  A(e)  have 


the  Smith  form  given  in  (3.2)  and  define 
y(t)  =  P_1(e)x(t) 


so  that 


Y(t)  =  D(e)A(e)y(t) 


(5.3) 


(5.4) 


where  A(e)  =  Q(e)P(e).  In  Section  3  and  4  we  performed  a  slightly  different 
similarity  transformation  and  also  replaced  A(e)  by  A  =  A(0) .  In  the  present 
context  we  cannot  throw  away  the  e-dependent  terms  in  A(e).  However,  as  the 
following  result  shows,  we  can  do  so  in  the  similarity  transformation  relating 
x(t)  and  y(t). 


Lemma  5.1:  Suppose  that  (A^(e) . A^e);  T)  defines  a  weak  time 

scale  decomposition  of  D(e)A(e).  Then  (A^(e) . An(e);  TP  *(0)) 

defines  one  for  A(e). 


Proof :  See  Appendix  D 

Let  us  introduce  some  notation.  Specifically,  let 


Aj  (e)  =  D(e)A(e)  = 


AnV)  R^fe)  ' 

eS  °(e)  tF  °(e) 


(5.5) 


where,  as  in  (4.3),  the  partition  indicated  is  compatible  with  that  of 

D(e)  =  diag( I ,  el . en_1I)  =  diag(I,  eD2(e))  (5.6) 

By  assumption  A^(e)  has  MSSNS,  so  A^^(e)  is  unimodular.  Consider  next  an 


arbitrary  (possibly  £-dependent)  matrix 


and  define  two  similarity  transformations  on  F: 


0F0  1  = 

I 

F  *F 

11  r  12 

F 

• 

I 

-F  ^F 

rll  12 

0 

I 

0 

I 

■I 

I 

0 

I 

0 

*F#  = 

-F  F  -1 

I 

F 

F  F  —  ^ 
Lr21r  11 

I 

L  21  11 

(5-7) 


(5.8) 


(5.9) 


We  also  define  a  thrid  similarity  transformation,  T  obtained  by  first  applying 
the  9-transformation  to  F  and  then  applying  the  'P  transformation  to  0F9  * 

(i.e.  we  construct  $  using  the  blocks  of  0F0  *).  We  can  now  state  the 
following. 

Lemma  5.2-  Define  the  following  sequences  of  matrices: 

r^i(fc)  =  T-transformation  for  A^(e)  (5.10) 

A1i+1(e)  =  r^eJA^fejr^e)'1  (5.11) 

where  Aj^(e)  is  given  in  (5.5).  Then  A^(e)  has  the  form 


A1i(e)  = 


A11  <£' 

i+U  i, 


£  Rx  (e) 


(5.12) 


e1  ^'(e)  eF^fe) 


where  A^fe)  is  unimodular.  Furthermore,  A^fe)  and  eF^e)  converge  to  the 
matrices  appearing  in  the  block-diagonal izat ion  of  A^(e)  obtained  as  in 
(4.21)  -  (4.26). 

Proof •  Equation  (5.12)  can  be  verified  by  direct  calculation.  See  [6,10,13] 
for  the  convergence  result  (which  is  not  used  in  what  follows). 

In  Section  4,  we  contented  ourselves  both  with  replacing  A(e)  in  (5.5) 
with  A(0)  =  A  and  with  performing  only  the  first  step  of  the  iteration.  In 
the  present  context  we  can  do  neither  of  these.  On  the  other  hand,  it  is 
still  not  necessary  to  go  to  the  limit.  To  make  this  precise,  we  begin  with 
some  notation.  Specifically,  let  N  denote  the  dimension  of  A(e);  X^e))  the 
eigenvalues  of  A(&);  and  M  an  upper  bound  on  the  maximum  order  of  the  real 


parts  of  the  X^Afe)).  i.e. 


0(Re[X.(A(e))])  <  M,  U1.....N 


(5.13) 


Since  we  have  assumed  that  A(e)  is  Hurwitz  such  a  bound  can  be  found.  For 
example,  if  A(e)  is  a  polynomial  matrix,  we  can  take  M  equal  to  the 
highest-order  power  of  e  appearing  in  |A(e)|. 


Given  N  and  M,  let 


K  =  NM  +  1 


(5.14) 


and  consider  carrying  out  K  steps  of  the  iteration  described  in  Lemma  5.2. 
This  produces 


UU.*J 
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A  K(  A 

A:  (e) 


Ai i <fe) 

K+lc  K+l 

6.  o , 


feKR1K(t) 

£A2°(fc) 


~  K  0  K 

where  we  have  introduced  the  notation  A^(e)  =  A^  (a)  and  A^  (a)  =  (a). 

Next,  we  perform  the  same  procedure  at  the  next  time  scale.  That  is,  write 


U) 


A22°(fc)  R2°(t) 
eS2°(e)  cF2°(£) 


(5.16) 


and  perform  K  steps  of  the  interaction  in  Lemma  5.2  involving  the  sequence 
Tg * ( e )  and  producing  A22(e)  and  A^(£).  Continuing  this  process  we  obtain  a 

<v  <v» 

complete  sequence  A^(e) . Ann^6^  can  state  t^le  following 

Theorem  5. 2-  Suppose  that  A(fc)  is  Hurwitz  and  has  MSSIES.  Then 
D(e)A(e)  has  a  weak  time  scale  decomposition  as  in  (5.2)  with  A.(e) 

=  A^e)  and  T  as  in  (4.17)  -  (4.19). 


Proof •  A  straightforward  calculation  shows  that 
2(£)D(£)A(£)2_1(£)  =  C(e)  +  H(e) 

where 

C(e)  =  diag(An(e),  eA99(e) . an~l\U))  (5.17) 


.  *'4.*4  .‘*.*4  «-<  t<.»  t.t  ijk  1 


0(H(e) )  =  K.  and 


2(e)  =  2jj(e)...2°(e)2n^(e)...2°(e)2^(e)...2°(e) 


(5. 18a) 


where 


£(‘>  -  r!fc) 


(5. 18b) 


k  1  0 

2  .(e)  =  , 

1  0  Tk(e) 

iv  ’ 


(5. 18c) 


As  in  Lemma  5.1,  we  can  replace  2(e)  by  2(0).  However  2(0)  =  T,  since 


T.(e)  =  I  for  k  >  0  and 
l  ’ 


I  A  °(0)-1R?(0) 

q(e)  =  11  1 

0  I 


(5.19) 


0  0  ~ 

with  A^  (0)  and  R.  (0)  equal  to  A^  and  Ri ,  respectively,  from  (4.14),  (4.15) 


What  remains  to  be  shown,  then  is  that  G(e)  and  G(e)  +  H(e)  are 


asymptotically  equivalent.  This  is  done  in  Appendix  E. 


The  key  idea  behind  this  result  is  that  we  must  approximate  the 


eigenstructure  of  A(e)  accurately  up  to  at  least  the  order  of  the  damping  in 


each  eigenmode.  For  example,  the  matrix 


-1  -fc-e 


Here  0(H(t))  denotes  the  minimum  order  of  all  i-lcments  of  H(e).  As  an  aside, 
note  that  the  diagonal  blocks  of  H(fc)  are  zero 


WVWVWVt 
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is  asymptotically  equivalent  to  the  matrix  in  (5.1)  —  i.e.  it  is  allowable  to 

2 

neglect  the  higher  order  (e  )  damping.  On  the  other  hand,  the  two  matrices 


r  2  , 

r  2  , 

-e  1 

-6  l-£ 

.-1  -e2 

* 

2 

L-l+fe  -e  J 

are  not  asymptotically  equivalent  since,  compared  to  the  order  of  damping,  the 
difference  in  frequency  (between  1  and  1+e)  is  very  significant. 

What  the  procedure  we  have  described  does  is  to  perform  a  sufficient 
number  of  iterations  to  guarantee  that  the  difference  between  the  eigenvalues 
of  A(e)  and  its  approximant  are  of  higher  order  than  the  real  (i.e.  the 
damping)  part.  Admittedly  the  procedure  is  conservative  —  typically  one  can 
get  by  with  fewer  iterations  and  can  discard  additional  higher-order  terms 
retained  by  the  procedure  —  but  it  is  guaranteed  to  work. 


6.  Assignment  of  Time  Scales  by  State  Feedback 

The  results  of  Sections  3  and  4  establish  the  role  of  the  invariant 
factors  of  A(e)  in  determining  the  time  scales  of  the  undriven  system  (1.1). 
For  the  driven  system  (1.11),  it  is  then  natural  to  pose  the  question  of  time 
scale  or  invariant  factor  assignment.  Specifically  it  is  of  interest  to 
determine  what  freedom  there  is  in  assigning  the  invariant  factors  of 

x(t)  =  F(e)x( t) ,  F(e)  =  A(e)  +  B(e)K(e)  (6.1) 

by  application  of  state  feedback  as  in  (1.12).  The  following  is  a  result  in 
this  direction. 

Theorem  6.1:  Assume  that  A(e),  B(e)  are  left  coprime,  i.e.  that 
[A(0)  B(0)]  has  full  row  rank.  Let  b  denote  the  rank  of  B(0). 

Then 

1.  F(e)  can  have  no  more  than  b  non-uni t  invariant  factors. 

jl  jb 

2.  There  exists  a  K(e)  such  that  F(e)  has  e  ,...t  as  its 
invariant  factors,  for  arbitrary  non-negative  integers 

00 

jj . (with  the  convention  that  e  =0). 

Proof :  We  first  show  that  we  can  further  assume  that 

A(t)  =  diag  (1 . 1.  e  1 . e  L)  ,  k.  >  0  (6.2) 

and  that  B(e)  is  upper  triangular.  Specifically,  suppose  that  A(e)  has  the 
Smith  form  given  in  (3.2).  We  can  then  write 

F(e)  =  P(e)[D(e)  +  P'1(e)B(£)K(e)Q"1(e)]Q(e)  (6.3) 

Thus  we  can  equivalently  consider  the  invariant  factors  of  D(e)  +  B(e)K(e), 
where  B(e)  =  P  *(e)B(e),  K(e)  =  K(e)Q  *(e).  Furthermore,  using  elementary 

_  A 

column  operations  we  can  show  that  B(e)U(e)  =  B(e)  where  U(e)  is  unimodular 
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and  B(e)  is  upper  triangular.  Consequently  we  have  the  equivalent  problem  of 

A  /"V  A  _  2  _ 

invariant  factor  assignment  for  D(e)  +  B(e)K(e),  where  K(e)  =  U  (e)K(e). 

Suppose  then  that  A(e)  is  given  by  (6.2)  and  B(e)  is  upper  triangular. 

Furthermore,  for  notational  simplicity  and  without  loss  of  generality  we 

assume  that  both  A(e)  and  B(e)  are  NxN.  Let  us  first  prove  the  second  part  of 

the  theorem  statement.  Note  first  that  for  [A(0) ,  B(0)]  to  have  full  row  rank 

9 

it  must  be  ture  that  L  <  b,  and  B  must  have  the  form 


B(e)  = 


L 


(6.4) 


where  *  represents  an  aribtrary  element  in  W. 

Assume  first  that  L  =  b.  Than  we  can  construct  a  unimodular  matrix  V(e) 
so  that 


f*  x  x 


B(e)V(e) 


}b 


(6.5) 


and  let 


K(e)  =  V(e)diag  (0 . 0.e 


-  6 


4  kb, 

,£  -  e  ) 


(6.6) 


^Actually  what  we  can  conclude  is  that  the  last  L  diagonal  elements  of  B(0) 
are  nonzero.  By  right-multiplications  we  can  make  these  values  unity. 
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It  is  straightforward  then  to  show  that  A(e)  +  B(e)K(fc)  has  the  desired 
invariant  factors.  If  L  <  b,  we  are  in  essence  replacing  some  of  the  unit 
invariant  factors  of  A(e)  with  nonunit  invariant  factors.  Since  rank  B(0)  = 
b.  b  -  L  of  the  first  N  -  L  columns  of  (6.4)  are  linearly  independent  at  e  = 
0.  Then,  just  as  in  constructing  (6.5),  we  can  construct  a  unimodular  matrix 
V(e)  so  that 


B(e)V(e) 


*  w . w 

0 _ 0  1  0 . 0 

0.  .0  1  0 . 0 

*  * . * 


0  I 


(6.7) 


i.e.  so  that  b-L  of  the  first  N-L  rows  are  zero  except  for  a  single  entry  of 
unity,  and  so  that  these  rows  and  the  last  L  rows  are  linearly  independent. 

In  this  case,  it  is  then  simply  a  matter  of  performing  a  permutation 
similarity  transformation  so  that  the  transformed  versions  of  A(e)  is  as  in 
(6.2)  with  some  of  the  =  0,  while  the  transformed  version  of  B(e)V(e)  is 
given  by  (6.5).  From  this  point  on  the  construction  is  the  same  as  before. 

To  prove  the  first  statement  in  the  theorem,  let  M  =  rank(A(0)  +B(0)K(0)) 
=  number  of  unit  invariant  factors  of  F(e).  Also,  assume  that  V(e)  has  been 
constructed  so  that  (6.5)  holds  (perhaps  after  the  permutation  similarity 
transformation  described  previously  if  L  <  b).  Letting  K(e)  =  V  ^(e)K(e),  we 
see  that 


A(0)  +  B(0)K(0)  = 
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where  the  *  terms  may  be  zero  or  nonzero;  however,  since  b  >  L,  there  are  b-L 
independnet  column  vectors  in  the  first  n-L  columns  of  the  matrix  multiplying 
K(0)  in  (6.8).  Consequently,  adding  B(0)K(0)  to  A(0)  can  reduce  the  rank  of 
A(0)  by  at  most  b-L.  Thus 

M  >  N  -  L  -  (b-L)  =  N-b  (6.9) 

Some  results  are  also  available  for  the  case  of  non-coprime  A(e),B(e). 

In  this  case  F(e)  is  of  the  form 

F(e)  =  W(t)F(e)  (6.10) 

where 

F(e)  =  A(e)  +  B(e)K(e)  (6-11) 

Here  W(e)  is  a  greatest  common  left  division  of  A(e).B(e),  and  A(e),B(e)  are 
left  coprime.  If  the  invariant  factors  of  F(e),  W(e),  and  F(e)  are  denoted  by 
fj(e.),  wi(e),  and  fj(e)  and  ordered  such  that  the  i1"*1  one  divides  the  (i+1)1"  , 
we  have  (thanks  to  the  Binet-Cauchy  formula  [14]) 

wi(a)|fi(t)  and  ^(OUjU)  (6.12) 

The  first  divisibility  condition  in  (6.12)  shows  that  every  invariant  factor 
of  F(e)  must  contain  the  corresponding  invariant  factor  of  W(e).  The  f.(e) 
are  governed  by  Theorem  6.1,  and  conclusions  about  the  f^(e)  can  then  be  drawn 
from  the  second  divisibility  condition  in  (6.12). 


7.  Conclusions 


In  this  paper  we  have  developed  an  algebraic  approach  to  time  scale 
analysis  of  singularly  perturbed  linear  systems  that  exposes  the  role  played 
by  the  Smith  form  of  A(e)  viewed  as  a  matrix  over  the  ring  of  functions 
analytic  at  e  =  0.  This  approach  bridges  the  gap  between  previous  easily 
interpreted  but  restricted  results  [1],  [2]  and  more  recent  results  [3]  that 
are  completely  general  but  quite  intricate.  Our  work  not  only  provides  a 
simple  interpretation  of  the  MSSNS  condition  introduced  in  [3]  in  terms  of  the 
inver tibi 1 i ty  of  successive  Schur  complements  of  a  particular  matrix  but  also 
allows  us  to  state  and  prove  a  strengthened  and  more  precise  version  of  the 
main  result  of  [3]  using  the  new  concept  of  a  strong  multiple  time  scale 
decomposi tion. 

The  framework  and  concepts  introduced  in  this  paper  also  open  the  way  for 
the  investigation  of  additional  questions.  Several  of  these  we  have 
considered  here  as  well.  In  particular,  we  have  investigated  the  relaxing  of 
the  so-called  MSST  condition  by  developing  a  procedure  involving  iterated 
Schur  complementation  in  order  to  guarantee  that  weak  but  essential  damping 
terms  are  retained.  In  addition,  we  have  investigated  the  problem  of 
time-scale  modification  via  state  feedback,  which  in  our  context  corresponds 
to  changing  the  invariant  factors  of  the  system  matrix.  Another  question  that 
can  be  asked  concerns  the  fact  that  the  Smith  decomposition  is  not  unique.  As 
shown  in  [8],  while  the  use  of  different  Smith  decompositions  leads  to 
different  time  scale  approximations,  the  successive  Schur  complements  in  these 
approximations  are  similar.  Also,  there  is  the  problem  of  computing  the  Smith 
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decomposition  of  A(e).  Some  ideas  related  to  this  are  given  in  [8],  but  these 
remain  to  be  developed.  In  a  sense  we  have  traded  the  difficult  tasks  of 
computing  e-dependent  projections  and  pseudo-inverses  that  are  needed  in  the 
approach  in  [3]  for  the  Smith  form  computation  in  our  approach.  However,  in 
our  work  this  computation  is  identified  as  a  separate  task  which  need  not  be 
carried  through  the  remaining  analysis  and  therefore  does  not  obscure  the 
intuition  behind  our  results. 

Finally,  note  that  in  [3]  the  orders  of  the  various  time  scales  of  (1.1) 
are  shown  to  correspond  to  the  orders  of  the  eigenvalues  of  A(e).  On  the 
other  hand,  in  this  paper  we  have  shown  that  the  orders  of  the  invariant 
factors  determine  the  time  scales.  It  should  not  come  as  too  much  of  a 
surprise  that  there  is  a  relationship  between  the  orders  of  eigenvalues  and 
invariant  factors  and  that  the  MSSNS  condition  plays  a  central  role  in  this 
relationship.  This  is  the  subject  of  a  forthcoming  paper. 
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Appendix  A :  Proof  of  Proposition  2.1 

This  proof  of  (2.17)  uses  several  of  the  ideas  introduced  and  developed  in 
Sections  3  and  4.  We  first  need  the  following 

Lemma  A. 1 •  Let  H(e)  be  obtained  form  A(t)  by  a  similarity 
transformation 


H(e)  =  S(fe)A(e)S'1(e) 
where  S(e)  is  unimodular.  Then 


(A- 1 ) 


(A. 2) 


This  result  follows  easily  from  the  fact  that  the  sequence  of  eigenprojections 
and  successive  system  matrices  defined  as  in  (2.12)  -  (2.16)  for  A(e)  and  H(e) 
are  all  related  by  the  same  similarity  transformation.  Equation  (A. 2)  then 
follows  on  examination  of  the  leading-order  terms  of  the  successive  system 
matrices . 

Consider  next  any  G(t)  fe  ?{A(e)},  i.e. 

G(e)  =  U(e)A(t)V(e)  (A. 3) 

wi  th 

U(0)  =  V(0)  =  I  (A. 4) 

Then  by  performing  similarity  transformations,  it  is  straightforward  to  check 
that  Proposition  2.1  will  be  proved  if  we  can  verify 

Lemma  A. 2-  Proposition  2.1  holds  if 

A(e)  =  D(e)A  (A. 5) 

C<£)  =  D(e)A(e)  (A. 6) 

with  A(0)  =  A,  which  is  invertible. 
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The  proof  of  this  lemma  is  a  straightforward  variation  on  the  development  j 

l 

in  Section  4.  As  in  Section  3,  let  us  assume,  without  loss  of  generality  that  i 

kl 

=  0  (since  otherwise  we  can  divide  (A. 5),  (A. 6)  by  e  ).  The  result  is 
then  proved  by  induction  on  n,  the  number  of  time  scales.  For  n=l  the  result 
is  immediate,  since 


A(e)  =  diag( I ,  OJA  = 


(A. 7) 


G(e)  =  diag( I ,  O)A(e) 


’All(t)  A12(£) 
0  0 


(A. 8) 


Clearly 


A11  A12 
.  0  0  . 


(A. 9) 


A  G 

Furthermore,  thanks  to  Lemma  4.1,  Fqq  =  Fqq  has  SSNS  if  and  only  if  A^  is 
invertible.  If  A^  is  singular,  the  procedure  stops.  If  A^  is  invertible, 
we  have  already  achieved  the  normal  rank  of  A(e)  (and  G(e))  so  that  all 
subsequent  F^q’s  are  equal  to  0.  In  either  case  the  lemma  is  verified. 

I f  n  >  1 ,  then 


and  G(e)  has  a  form  analogous  to  (A. 10)  wiht  A.  .(e)  replacing  A...  Again 

(A. 9)  holds,  and,  as  before,  the  procedure  stops  if  A^  is  singular.  If  A^ 

is  invertible,  we  can  use  the  same  procedure  as  sketched  at  the  end  of 

Section  4  to  block  diagonalize  A(e)  and  G(e).  Specifically,  consider 

equations  (4.21)  and  (4.22),  where,  we  replace  A^  by  A^,  by  A^2>  Sj(e)  by 

DM(e)A2j.  and  F^e)  by  D^(e)A22.  Again  because  of  the  invertibility  of  A^, 

solutions  LA(e)  and  H^(e)  exist  to  these  equations,  with 

L1A(0)  =  -  An_1A12  HjV)  =  DM(0)A21An_1  (A.  12) 

Similarly,  we  can  solve  (4.21),  (4.22)  with  analogous  replacements  but  with 

G  G 

A_(e)  substituted  for  A^.  This  yields  solutions  L ^  (e),  (e).  Applying 

the  corresponding  diagonalizing  similarity  transformations  (4.24),  (4.25)  to 

A  G 

A(e)  and  G(e),  and  noting  that  (0)  =  (0),  we  see  that,  thanks  to  Lemma 

A.l,  we  have  reduced  the  problem  to  one  with  one  fewer  time  scale  —  i.e.  we 
are  lef  t  to  examine 
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Appendix  B:  Proof  of  Proposition  2.2 

Without  loss  of  generality  we  assume  that  the  similarity  transformation  T 

in  (2.1)  is  the  identity  —  if  this  is  not  the  case  we  can  simply  perform  an 

initial  e-independent  similarity  transformation  on  Aq(c).  Furthermore,  since 

A, . A  are  assumed  to  be  semistable,  we  can  perform  another  e-independent 

1  n 

similarity  transformation  so  that  what  we  are  given  are  Hurwitz  matrices 


G. . G  so  that 

1  n 

An(e)t  k1  k 

lim  sup||e  -exp(diag  (0,  e  G. . e  nG  )t}||  =  0 

eiO  t>0 


(B.l) 


and  what  we  would  like  to  show  is  that 


diag  (0 . 0,  G.,0 . 0)  k=k^ 


(B  -  2 ) 


otherwise 


As  a  first  step,  note  that  (B.l)  implies  that  for  any  integer  r 


An(e)t  k.  k. 

lim  sup  _  |  |e  -exp(diag  (0,  e  G1 . e  JG  ,0 . 0}  1 1  ]  =  0  (B.3) 

elO  te[0,e  r)  J 


where 


kj  <  r  <  kj+1  ,  j=0 . n 


(B.4) 


(here,  for  completeness  kQ  =  0,  kn+j  =  ")•  Note  also  that,  since  F00  ~  V°> 


An(e)t  Fmt 

lim  sup  | | e  -  e  |  |  =  0 

eiO  te[0, 1) 


(B.5) 
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From  (B.3)  -  (B.5)  we  can  conclude  that  if  >  0,  Fqq  =  0,  Pq(£)  =  I.  and 
A^(t)  =  Aq(£)/£.  Consequently,  we  can  simply  replace  A^e)  in  (A.l)  with 
Aj(e)  and  reduce  each  of  the  by  1.  Continuing  in  this  fashion  we  find  that 

Fj^q  =  0,  k  <  kj .  From  [3]  we  then  have 

An(fc)t  k 

lim  sup  | | e  -  exp{F,  „£  t} j |  =  0  (B.6) 

elO  -k  1 

tfe[0,£  l) 

and  from  (B.3),  (B.6)  we  conclude  that 

Fk  =  diag(0.  Gr0 . 0)  (B.7) 

The  remainder  of  the  proof  proceeds  by  induction  on  n.  The  case  of  n=l 
is  essentially  complete,  since  in  this  case  the  sup  on  the  left-hand  side  of 

p 

(A. 6)  can  be  taken  over  [0,  1/e  )  for  any  r  >  k^ .  Consequently  an  argument 
identical  to  the  one  used  in  the  preceding  paragraph  shows  that  F^  =  0  for 
all  k  >  kj .  To  consider  the  case  of  n  >  1,  we  assume,  without  loss  of 
generality,  that  k^  =  0  (since  as  we  have  seen,  if  k^  >  0  then  Aq(e)  is 
kl 

divisible  by  e  so  we  can  rescale  time  to  eliminate  this  factor).  Next, 
write  Aq(c.)  as  the  sum  of  two  commuting  matrices 
A0(O  =  P0(e)A0(e)  +  [I-P0(e)]A0(6) 

=  eA^e)  +  [I  -  P0(e)]A0(e)  (B.8) 

Note  that,  from  [3]  and  (B.7) 

[I-P0(e)]A0(e)t 

lim  sup||e  -  exp(diag(0,  G. ,  0 . 0) t ) | |  =  0  (B.9) 

elO  t>0 
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Then,  using  (B.8)  and  performing  several  standard  manipulations  we  obtain  the 
following 


A0(e)t 


k9  kn 

2^-,  rv 


exp{diag  (O.G^.e  G2 . e  G^)t} 


Aj(fc)fct 


<  e 


k9  k_ 

2.^  n^ 


-  exp{diag(0,0,e  G2 . e  Gn ) t } j  | .  |  | e 


[I-P0(e)]A0(e)t 


[I-P  (e)]A  (e)t 

+  I  |e  -  exp{diag(0,G1 ,0 . 0)t}||. 


k  k 

| |exp(diag(0,0,e  2G2 . e  nGn)t}||  (B.10) 


Note  that  since  n  >  1,  (B.9)  implies  that 


[I~P0(fc)]A0(e)t 


is  bounded  away  from  zero  uniformly  in  t.  Consequently  (B.l),  (B.9),  (B.10), 
and  the  semistability  of  Gg . G^  imply  that 


k  -1 


A  (e)t  k„-l 

lim  sup | |e  -exp(diag  (0,0, e  G„ . e  n  G)t}||=0 

e  10  t  >0  Z  n 


(B. 11) 


and  consequently  (B.2)  follows  by  induction. 
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Appendix  C:  Verification  of  Equation  (4.13) 

Let  us  rewrite  (4.11)  as 


C(e) 


Gll(£'  G12(e) 

G21 (e)  G22(6). 


where 


Gll(e) 

gi2(-) 

G2i(-) 

G22(e) 


An  +  e  2An  1R1D2(e)W1 


(C.l) 


(C. 2a) 
(C. 2b) 
(C.2c) 
(C.2d) 


Note  that  G^e)  is  invertible  in  a  neighborhood  of  e  =  0.  Let 


C^e) 


I 

^21 C & )Gi i  (e) 


0  ' 
I  . 


I 

i 

I  E(e) 


I 

0 


-Gjl  \e)G12{e) 
I 


(C.3) 


(C.4) 


From  (C.2)  we  see  that  C^(0)  =  E(0)  =  I,  and  a  straightforward  calculation 
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Appendix  D:  Proof  of  Leaaa  5.1 

We  have  that 

lim  sup  ||eD(e^&)t  -T-1exp{diag[An(e)  .eA.  (e) . en  lk  ,  (e)]t}T||  =  0 

elO  t>0 

(D.l) 

Therefore 

lim  sup | |eA^6^ t-P(e)T  1exp{diag[A_(e) ,eA1 (e) . en  *A  ,(e)]t}TP  1 (fc)  |  |  =  0 

elO  t>0 

(D.  2) 

What  we  must  show  is 

lim  sup  |  |P(e)T  *exp{diag[An(e) . e.n  ^A  ..(fc)]t}TP  *(e) 

elO  t>0  n-A 

-  P(0)T-1exp{diag[A0(e) . e^A^eJjtyr^O))  |  |  =  0  (D.3) 

A  simple  triangle  inequality  argument  shows  that  the  left-hand  side  of  (D.3) 
is  bounded  above  by 

lim  sup|  |(P(e)-P(0))T  ^xpfdiagfAfTe) . en  *A  (e)]t}TP  1  (e.)  |  |  + 

eiO  t>0 

lim  sup  | |P(0)T_1exp{diag[A  (e) . e^A. (e)]t}T(P_1(£)-P_1(0)) | |  (D.4) 

£ 10  t>0 

The  first  term  in  (D.4)  is  in  turn  bounded  above  by 

lim{  |  |P(£)-P(0)  1  |  |  |P_1(e)  |  M|T|  M  |T-1 1  I 
£  10 

•  sup| |exp(diag[A  (£) . £n  *A ,(e)]t}||  (D.5) 

t>0  u  n  1 

From  the  construction  in  Section  5,  we  know  that  each  Aj(e)  is  Hurwitz 
for  £  >  0  and,  since  A^(0)  has  MSSIES,  we  know  that 

| |exp{diag[AQ(0) . en  *An_^(0)]t} | |  is  bounded.  From  this  we  conclude  that 
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Appendix  E:  Completion  of  the  Proof  of  Theorem  5.2 

The  result  we  need  to  prove  is  the  following 
Proposition  E. 1 :  Suppose  that  the  NxN  matrix  G(e)  is  Hurwitz.  Suppose 
further  that 

0[Re(X.(G(e))]  <  M  i=l . N  (E.l) 

and  let  K  =  MN  +  1.  Then  G(e)  is  asymptotically  equivalent  to  G(e)  +  H(e), 
where  H(e)  is  any  matrix  with  0(H(e))  =  K. 

Proof :  The  proof  is  a  variation  on  the  methods  in  [3.  7,  8].  First  from  [7] 
we  have  the  following 


Lemma  E.l:  Let  A(e)  =  B(e)  +  £:  C(e)  be  an  NxN  matrix.  Then 
min  0[A . (A(e) )  -  A.(B(e))]  >  p/N  (E.2) 

.1  l 


*  i . 


(E.3) 


Consequently  in  our  case 

0[A.(G(e))  -  X.(G(e)  +  H(e))]  >  0[Re(X . (G(e) ) ) ] 

Next,  recall  the  definition  of  the  resolvent  of  a  matrix  A(e) 

R(X.A)  =  [A(e)  -  XI]"1  (E .4) 

so  that 

A(e] 

2iri 

k  ‘  k 


(E.5) 


where  the  are  positively-oriented  contours  enclosing  disjoint  portions  of 

the  complex  plane  and  all  of  the  eigenvalues  of  A(e).  Consider,  then 

2Tri[e^G^e)+H^£^t  -  eG{6^]  =  ^  |  eXt[R(X,G)  -  R(X,G+H)]dX  (E.6) 

k  fk 

where  we  choose  the  carefully.  Specifically  is  a  circle  centered  at 


*  •  * ■  v  1 4  m\  % v  r*/V  v  v  v  v.V-V*  A* ^ /a’** , «/■  Jfii 
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X^(G(e)),  of  radius  of  order  0[Re(X^(G(e)))] ,  and  completely  contained  in  the 
left-half  plane  {Re(A)  <  0}.  More  precisely,  we  require  the  maximum  value  of 
Re(A)  on  to  also  be  of  order  0[Re  (X^ (G(e ) ) ) ] . ^  Also,  for  e  small  enough 
(E.3)  guarantees  that  this  circle  includes  X^(G(e)  +  H(e)).  The  circle  may 
also  include  other  pairs  of  eigenvalues,  but  for  e  sufficiently  small  this 
happens  only  if 

0[X1(G(e))-X.(G(e))]  <  min  {0[Re(X1 (G(e ) ) ) ] .  0[Re(X  .  (G(e) ) ) ] }  (E.7) 

Consider  next  a  single  term  in  (E.6)  and  suppose  that  the  radius  of  is 

of  order  m.  If  we  let  A‘  =  A/e™  we  can  rewrite  this  term  as 
nu  , 

e£  t[R(emA '  ,G)  -  RfeV,  G+H) ]emdA 1  (E.8) 

where  ,  the  image  of  under  this  mapping,  has  radius  of  order  1,  is 
completely  contained  in  the  left-half  plane,  and  in  fact  consists  of  points 
with  negative  real  parts  of  order  1.  Consequently,  the  norm  of  (E.S)  is 
bounded  above  by 

llRfeV.G)  -  R(eV  ,G+H)  |  |emdA'  (E.9) 

k 

Also,  we  can  write 

R(A,G)  -  R(A,G+H)  =  R(X,G){I  -  [I+HR(X.G)]_1 }  (E.10) 

m  k 

Note  that,  thanks  to  (E.3)  and  (E.7),  R(X,G)  is  of  order  1/t  on  V  . 
Consequently  (since  m  <  M)  HR(X,G)  is  of  order  at  least  m(N-l)  +  1.  and  we  can 
write  the  series 


^For  example,  the  circle  (X: |X-X^(G(t)) |  =  Re(X  (G(e)))}  will  due  unless 
another  singularity  lies  on  it. 


(E.l 


R(A,G)  -  R(A,G+H)  =  R(A,G)  J  (-l)n  [HR(A,G)]n 

n=l 

which  converges  uniformly  for  A  £  Obviously  the  same  statements  cam  be 

made  for  R(emA'  ,  G)  and  HRfe’V,  G)  on  iy  .  and  therefore  we  conclude  that 
0(  |  |R(emA‘  ,G)  -  R(e.nVv  '  ,G-»H)  |  |em)  >  m(N-l)  +  1  (E.l 

uniformly  on  iy  .  Since  iy  has  perimeter  of  order  1  in  length,  (E.9) 
converges  to  0  as  elO,  and  the  result  follows. 


